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Introduction 

A major challenge in the quest to constrain the accretion his- 
tory of the Universe is to compile unbiased samples of active 
galactic nuclei (AGN). Obscuration in the vicinity of the super- 
massive black hole (SBH) and dilution of the AGN emission 
from stellar light associated with the host galaxy (e.g. Moran 
et al. 2002, Severgnini et al. 2003), introduce severe selection 
effects against certain types of AGN. X-ray observations, espe- 
cially at energies above about 2 keV, have been shown to be ef- 
ficient in minimising these biases. Hard X-rays, 2-lOkeV, can 
penetrate high column densities of gas and dust clouds and suf- 
fer minimal contamination by emission from the host galaxy, e.g. 
X-ray binaries and hot gas (Hornschemeier et al. 2003). 

Recent observations at X-rays have highlighted the effec- 
tiveness of this wavelength regime in locating accreting SBHs. 
Firstly, X-ray luminous AGN have been revealed at the centres 
of apparently normal galaxies, with optical spectra dominated 
by stellar emission (e.g. Comastri et al. 2002; Georgantopoulos 
et al. 2003; Georgantopoulos & Georgakakis 2005). Secondly, 
the deepest X-ray survey to date, the 2 Ms Chandra Deep Field 
North (CDF-N), finds an AGN density of over ~ SOOOdeg ^ 
(Alexander et al. 2003; Bauer et al. 2004), an order of magnitude 
higher than optically selected broad-line QSOs, which reach 
number densities of only 400 deg~^ at /? = 24 mag (IWolf et al.L 

Eooj. 

The X-ray surveys in the 2-lOkeV band , although efli- 
cient in compiling nearly unbiased AGN samples are not per- 
fect either In the case of very high obscuring gas column den- 
sities, > lO^'^cm"^ , Compton scattering by bound electrons 
becomes important. This mechanism blocks photons with en- 
ergies below almost lOkeV at the rest frame. If AGN lie be- 
hind such high columns densities, one can only hope to de- 
tect them indirectly in the 2-lOkeV band. For example in the 
case of a toroidal geometry of the obscuring material X-ray 
photons are expected to be reflected off the inner surface the 
torus into our line of sight. This reflected component how- 
ever, is typically about 2 dex fainter than the direct (obscured) 
power-law emission. This causes even the 2-lOkeV surveys to 
under-sample the Compton-thick AGN population (To zzi et al., 
l2006t [Georgantopoulos. Georgakakis & Akylasl l2007h . A large 
number of active SBH are believed to belong to this class 
of heavily obscured sources. Up to 50% of the Seyfert- 
2s at low redshift for example, are C ompton-thick candi- 
dates JMaiolino & Rieke, 19951 iRisaliti. Maiolino & Salvati, 
119991; iGuainazzi. Matt & Perolal l2005t iPanessa et all 120061k 
Moreover, the peak of the diffuse X-ray background (XRB) 
spectrum at 30-40 keV (Gruber et al. 1999, Frontera et al. 2007, 
Churazov et al. 2007) can only be modeled by invoking a popu- 
lation of Compton-thick AGN, which are missing from existing 
2-IOkeV X-ray surveys (Gilli et al. 2007). The evidence above 
underlines the significance of Compton-thick AGN for under- 
standing both the XRB and the accretion on SBHs across the 
history of the Universe. The precise fraction of these heavily ob- 
scured sources is still uncertain however, by up to a factor of few 
(e.g. GilH et al. 2007). 

The bias of the 2-lOkeV X-ray imaging surveys against 
heavily obscured AGN has motivated alternative methods for 
finding such systems. The mid-infrared wavelength regime (mid- 
IR; 3 -30//m) in particular, has received much attention recently 
and is proposed as a potentially powerful tool for identifying 
deeply buried AGN. The optical and UV photons from the AGN 
heat the dust of the obscuring material around the central engine 
and are re-emitted as thermal radiation at the mid-IR. Therefore, 
surveys with the Spitzer infrared mission provide an excellent 



opportunity to recover these heavily obscured sources. The dis- 
advantage of mid-IR surveys is that the AGN are outnumbered 
by star-forming galaxies and some selection method is necessary 
to separate the two populations. It has been suggested that mid- 
IR colour-colour plots can be used to effectively isolate obscured 
AGN. The principle behind this approach is that luminous AGN 
have power-law spectral energy distribution (SED) in the mid- 
IR, while galaxies have characteristic black-body spectra that 
peak at about 1.6;um. The mid-IR colours of AGN are therefore 
redder than those of galaxies defining a characteristic wedge. 
Different co mbinations of mid -IR colours have been proposed to 
select AGN. 'Lacy et al.' (2004') and'Hatziminaoglou et"al1 (|2005|) 
define their selection window based on the mid-IR colours of lu- 
minous high-redshift QSOs identified in the Sloan Digital Sky 
Survey. These techniques were tested against control samples 
and turned to be particularly efficient when combined with an 
additional criterion. Lacy et al. (2004) impose a flux larger than 
1 mJy at 8.0 ^m. Stern et al. (2005) impose a magnitude limit 
of R=19 finding that their mid-IR colours selection technique 
provides a 90% efficiency in detectin g AGN. 

In a complementary approach Al onso-Herrero et al] (l2006l) 
and iDonlev et al.l (12007 5 select sources with power-law SEDs 
in the mid-IR. The X-ray identification rate of these sources 
is about 50% in t he 1 Ms Chandra Deep Field South 
(iGiacconi et al.L l2002h . increasing to a; 85% in the deeper 2 Ms 
CDF-N. Parallel to the mid-IR studies above, selection meth- 
ods using a combination of IR, optical and/or radio criteria have 
been developed, also claiming the identification of type-2 AGN 
(e.g. Martinez-Sansigre et al. 2005; Daddi et al. 2007; Fiore et 
al. 2008). 

The selection methods above have been highly influential by 
providing samples of deeply buried AGN candidates, some of 
which may be so obscured that remain undetected to the lim- 
its of the deepest current X-ray surveys. Even for these sources 
however, getting a handle on their X-ray properties is essential, 
firstly to control any residual contamination by non-AGN and 
secondly to estimate their contribution to the XRB. Stacking 
analysis is a powerful tool to study the mean X-ray properties 
of sources below the detection threshold of X-ray observations 
(e.g. Daddi et al. 2007; Fiore et al. 2008; Georgakakis et al. 
2008). Here we use X-ray stacking to assess the efficiency of 
mid-IR based selection methods in identifying obscured AGN. 
The data are from the northern field of the Great Observatories 
Origins Deep Survey (GOODS; Dickinson et al. 2003). We 
adopt Ho = 75 kms' Mpc', Dm = 0.3, Qa = 0.7. 

1 . The Data 

The CDF-N is centred at a = 12''36'"49'.4, 5 = -62°12'58" 
(J2000) and has been surveyed extensively over a range of wave- 
lengths by both ground-based facilities and space missions. The 
multiwaveband data in this field include Chandra X-ray obser- 
vations, Spitzer mid-IR photometry, HST/ACS high resolution 
optical imaging, deep optical photometry and spectroscopy from 
the largest ground based telescopes. 

The 2Ms Chandra survey of the CDF-N consists of 20 in- 
dividual ACIS-I (Advanced CCD Imaging Spectrometer) point- 
ings observed between 1999 and 2002. The combined observa- 
tions cover a total area of 447.8 arcmin^ and provide the deep- 
est X-ray sample c urrently available. Here we use the X-ray 
source catalogue of [Alexander et al.l (l2003l) . which consists of 
503 sources detected in at least one of the seven X-ray spectral 
bands defined by these authors in the range 0.3 - 10 keV. The flux 
limits in the 0.5 - 2 and the 2 - lOkeV bands are 2.5 x 10"'^ and 



1.4x10 erg cm s , respectively. The Galactic column den- 
sity towards the CDF-N is 1.6 x 10^'^ cm"^ dPickev & LockmaiJ, 
[1990). 

The most X-ray sensitive part of the CDF-N has been ob- 
served in the mid-lR by the Spitzer mission (IWerner et al.Ll2004b 
as part of the GOODS. These observations cover an area of 
about 10 X 16.5 arcmin^ in the CDF-N using both the IRAC 
(3.6, 4.5, 5.8 and 8.0 fim) and the MIPS (24 //m) instruments on- 
board Spitzer Here we use the 2nd data release (DR2) of the 
IRAC super-deep images (version 0.30) and the interim data re- 
lease (DRl-i-)of theMIPS 24 /im mosaic (version 0.36) provided 
by the GOODS team (Dickinson et al. 2003; Dickinson et al. 
in preparati on). Sources are d e tected in these images using the 
SExtractor (iBertin & Arnoutsl 1 19961) software. The complete- 
ness limits i.e. the flux where the counts peak before dropping 
off" because of incompleteness are 0.7, 0.7, 2.8, 3.5 and 56 yuJy 
in the 3.6, 4.5, 5.8, 8.0 and 24 jjm bands respectively. Full de- 
tails on the source extract ion and flux derivation are presented in 
iGeorgakakis et al.l (|2007[) . Much of the analysis presented here 
uses the 24yum selected sample, which consists of 1619 sources 
down to a (3cr) limiting flux of about 15 fi3y. 

Multi-waveband optical imaging (BVRIz) in the CDF-N re- 
gion has been obtained using the SUBARU 8.2-m telescope. 
U-band observations have been obtained at the 4-m Kitt Peak 
Telescope while HK' observations have been obtained at the 
2.2-m University of Hawaii telescope. All the observations are 
described in Capak et al. (2004). Deep HST ACS observa- 
tions in the B,V,I,Z bands are also available (Giavalisco et al. 
2004). However, we choose to use the catalogue of Capak et 
al. (2004) as this is clearly more suitable for the derivation of 
photometric resdhifts. Here, we use the 7?-band selected sample, 
which contains 4745 1 sources down to a limiting magnitude of 
Rab - 26.6 mag (5cr). These observations cover about 0.2 deg^ 
and extend beyond the GOODS field of view {k. 0.05 deg-). The 
24yum catalogue is first cross-correlated with the IRAC 3.6jjm 
catalogue using a matching radius of 2 arcsec. For 24yum sources 
with 3.6/im counterparts (95 per cent of the 24/im population), 
the more accurate 3.6^m positions are used to search for opti- 
cal identifications within a radius of 2 arcsec. There are 1409 
24yum/optical associations. Given the search radius of 2 arcsec 
and the distribution of the positional off'sets between the optical 
and the mid-IR sources, we expect a spurious identification rate 
of about 4.5% at the limit of the optical catalogue, /? a; 26.5 mag. 
A total of 336 out of the 1409 sources are excluded from the 
analysis, because they lie either close to a bright source, or at the 
edge of the surveyed area or they are not detected in at least one 
of the four IRAC bands. 

Optical spectroscopy in the CDF-N and in the GOODS- 
North is available from either observations that specifically tar- 
get the X-ray population in these fields (e.g. Barger et al. 2003, 
2005; Cowie et al. 2004) or the Keck Treasury Redshift Survey 
(TKRS). This program uses the DEIMOS spectrograph (Faber 
et al. 2003) at the Keck-II telescope to observe optically selected 
galaxies to Rab ~ 24.3 mag in the GOODS-North. The publicly 
available catalogue consists of about 1400 secure redshifts, of 
which about 500 are associated with 24/im sources. 



2. The sample selection 

The sample of mid-IR wedge AGN candidates is compiled from 
the GOODS-North 24yL(m catalogue. Two different AGN identifi- 
cation methods are emp loyed. The first one is based on the IRAC 
band colour criteria of lStern et al.l (l2005h . The second method 



combines optical and mid-IR wavebands by adopting criteria 
similar to those proposed by Fiore et al. (2008). 

2.1. The mid-IR wedge AGN candidates 

The mid-IR wedge AGN candidates lie in the area of the colour- 
colour diagram proposed by Stern et al. (2005lj. This selection 
method is graphically shown in Fig. [1] where in the upper left 
panel we plot the X-ray detected sample (173 sources) while in 
the upper right panel we plot all the sources without X-ray de- 
tection (896). Only lA/um sources with optical counterparts and 
detections in all 4 IRAC bands are employed. The sample of 
mid-IR wedge selected candidate AGN consists of 177 sources, 
46 of them with X-ray detection and 131 not individually de- 
tected i n the X-rays . 

The I Stern et al] (|2005|) criteria are preferred here over alter- 
native AGN identification methods which use mid-IR colours 
only, because they are likely to suffer less contamination from 
normal galaxies, while at the same time remain sensitive to ob- 
scured AGN (e.g. Barmby et al. 2006). 

For comparison, we also explore the X-ray properties of the 
sources outside the Stern colour wedge. This sample has 896 
24/vm sources with optical identifications and detections in all 4 
IRAC bands. 



2.2. Tine red/opticaiiy-faint mid-IR AGN candidates 

Fiore et al. (2008) proposed a method which combines optical, 
near-IR and mid-IR data to identify AGN. These authors used 
red optical colours (R - K > 4.5) in combination with high 
mid-IR (24jum) to optical (/?-band) flux ratio, /24^m//R > 1000, 
claiming a large number of Compton-thick AGN in the CDF- 
South. 

We use similar selection criteria to compile the sample of 
red/optically-faint mid-IR AGN candidates in the CDF-North. 
This consists of 103 24//m sources with 3.6;um detection, 
fiAiimlfR > 1000 and Rab - m^.e > 3.7 mag, where TO3.6 is the 
magnitude at 3.6/im in the AB system. The Fiore et al. (2008) 
R - K > 4.5 mag colour cut has been replaced by the nearly 
equivalent Rab - '«3.6 > 3.7 mag. This cutoff is adopted to select 
sources redder than the elliptical galaxy template of Coleman, 
Wu & Weedman (1980) a.tz~ 1. Optical counterparts are avail- 
able for 31 of these sources. Therefore most of our points refer 
to Rab - OT3.6 and fi^^mlfR lower limits. 

The position of the red/optically-faint mid-IR AGN candi- 
dates on the Stern et al. (2005) colour diagram is shown in Fig. 
[1] The majority of these sources lie outside and on the left of the 
AGN wedge on this diagram. 

3. The Spectral Energy Distribution fitting method 

For 24/im sources with optical counterparts (1409), we model 
the observed UV/optical to mid-IR Spectral Energy Distribution 
(SED). This is both to determine photometric redshifts for 
sources without optical spectroscopy and to get an estimate of 
the total infrared luminosity. 

The adopted template fitting method is that of Rowan- 
Robinson et al. (2005). This is a two step process. First, the 
L'^-band to 4.5 /im photometric data are fit using a library of 8 
templates described by Babbedge et al. (2004), six galaxies (E, 



' ([5.8] - [8.0]) > 0.6, ([3.6] - [4.5]) > 0.2([5.8] - [8.0]) + 0.18, and 
([3.6] - [4.5]) > 2.5([5.8] - [8.0]) - 3.5. [3.6], [4.5], [5.8], and [8.0] are 
the (Vega) magnitudes in the respective IRAC wavebands 



in 

I 0.5 

'to' 





1 ' ' 


' 1 




1 


1 1 1 




~ 




. 


• 1 






~ 








1 










• 


, 


/• ". 1 










• ■ . 




■'•• ' 










• a.. 




1 






" 


- 






••. 1 


• 




- 




t 




' J . 


•« 










'.' 


.^'''- 








- 


• 


\^ 


^^.' 






- 


_ 




•^ 


, 




* 


_ 


. 




• 






. 


_ 


•:;.*,?••. 




• 


•• 


_ 












_ 


. •.■* 








_ 


_ 


• 








_ 


"l 1 


1 1 1 1 1 1 


III 


1 


1 1 1 


1 1 1 




in 



I 0.5 



CO 




1 2 

[5.8] - [8.0] 



I 0.5 



1 2 

[5.8] - [8.0] 




1 2 

[5.8] - [8.0] 



Fig. 1. Mid-IR colour-c olour plot. For all c larity we use different panels. In all of them the wedge defines the AGN region following 
the selection criteria of Ste rn et al.l (|2005). Top left: the 24jt/m sources with X-ray counterparts; Top right: the 24yum population 
with detections in all 4 IRAC bands, not detected in X-ray wavelengths; bottom left: the red/optically-faint mid-IR AGN candidates 
(the crosses and the filled symbols represent the X-ray and non X-ray detections respectively); bottom right: colour tracks for 
different template SEDs as function of redshift up to z = 3. The dotted (blue) line represents a QSO template (Elvis et al. 1994), the 
continuous (magenta) Mrk273 (ULIRG, Seyfert 2), the dashed (green) Arp220 (ULIRG, starburst). For each template the position 
of the redshift z = is marked with a large dot. Small dots correspond to z=l steps. All templates are adapted from Donley et al. 
(2007) 



Sab, Sbc, Scd, Sdm and sb) and two AGN. One AGN template is 
based on the SDSS median composite quasar spectrum (Vanden 
Berk et al. 2001); this is extended to wavelengths longer than 
8555A using the IR part of the average QSO spectrum defined 
in Rowan-Robinson et al. (2004). In addition to the composite 
SDSS template, a simpler template has been used, based on the 
mean optical quasar spectrum of Rowan-Robinson (1995) span- 
ning the range 400A to 25 jum. Photometric redshifts are deter- 
mined at this stage. The accuracy is &/(l H- Zspec) ~ 0.04 (Icr 
rms), estimated based on the 513 sources with available optical 
spectroscopy. 

At the next step the mid-IR SED of the sources in the sam- 
ple is fit. At longer wavelengths (5.8 - 24jum) any dust may 
significantly contribute or even dominate the observed emission. 
Before fitting models to these wavelengths the stellar contribu- 
tion is subtracted from the photometric data by extrapolating the 
best-fit galaxy template from the previous step. The residuals 
are then fit with a mixture of four templates: cirrus (Efstathiou & 
Rowan-Robinson 2003), AGN dust tori (Rowan-Robinson 1995; 



Efstathiou & Rowan-Robinson 1995), M82 and Arp220 star- 
bursts (Efstathiou et al. 2000). The modeling above provides 
an estimate of the total infrared luminosity, Ltot, in the wave- 
length range 3 - 1000/im. As discussed by Rowan-Robinson et 
al. (2005) this is expected to be accurate within a factor of two. 

4. X-ray stacking 

4.1. Method 

Stacking techniques have been widely used in X-ray Astronomy 
to study the mean properties of source populations selected to 
have certain well defined properties a nd which are too X-ray 
faint to be detecte d individually (e.g . [Alexander et al.L 120011; 
iNandra et all 120021; lGeorgakakisetal.[l2003h . 

A fixed radius aperture is used to extract and to sum the X- 
ray photons at the positions of the various 24 fim selected sub- 
samples. Sources that lie close to or are associated with an X-ray 
detection are excluded from the analysis to avoid contamination 
of the stacked signal from the X-ray photons of detected sources. 
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Fig. 2. The redshift distribution of the mid-IR wedge selected 
AGN: X-ray detections (shaded) and non X-ray detections 
(open) histogram. 



Fig. 3. The IR luminosity distribution of the mid-lR wedge se- 
lected AGN: X-ray detections (shaded) and non X-ray detections 
(open) histogram. 



We adopt an extraction radius of 3 arcsec. This is found to max- 
imise the signal-to-noise ratio of the stacked signal. A 3 arcsec 
aperture encloses more than about 90% of the photons in the 
1 .5 - 4 keV spectral band at an off-axis angle of 5 arcmin. The 
average off-axis angle for our sources is 5.5 arcmin. 

The significance of the stacked signal depends on the value 
of the background. This is estimated using the smoothed back- 
ground maps produced by the WAVDETECT task of CIAO by 
summing the X-ray photons in regions around each source used 
in the stacking. The significance of the stacked signal in back- 
ground standard deviations is estimated by (T - B)/ Vb, where 
T and B are the total (source + background) and background 
counts respectively. 

An estimate of the mean spectral shape of the detected sig- 
nal is obtained by performing the stacking in two energy bands, 
soft (0.3 - 1.5 keV) and hard (1.5 - 4.0 keV). Fiore et al. (2008) 
has shown that this choice of bands minimises the instrumen- 
tal background contamination and maximises the sensitivity of 
the stacking analysis. Fluxes are determined by multiplying the 
stacked count rate by the appropriate energy conversion factor, 
which is estimated separately for each class of sources, based on 
the spectral shape of the stacked signal. 

4.2. Mid-IR wedge AGN 

Table [T] summarises the fraction of X-ray detected AGN in dif- 
ferent 24//m selected samples. The mid-IR wedge AGN candi- 
dates have a relatively high X-ray identification rate, about 26%. 
In contrast, only about 14% of the 24yL(m sources outside that 
wedge are associated with X-ray AGN. Fig. 2 shows that the 
mid-IR wedge AGN candidates lie at moderately high redshift 
on average, z ~ 1 .4. The X-ray detected members of this pop- 
ulation have a flatter redshift distribution that extends to z > 2. 
It appears that X-ray observations are more efficient in detecting 



Table 1. X-ray detected AGN fraction in different mid-IR sub- 
samples 



Sample 
(1) 



Ntot 
(2) 



Nx 
(3) 



A 
(%) 
(4) 



wedge mid-IR AGN 
24/im sources outside 
wedge 

red/optically-faint 
mid-IR AGN 



177 
896 

103 



46 
127 

20 



26 
14 

19 



The columns are: (1): sample definition. (2): NjoT is the total number of 24/j m sources. 
(3): N x is the number of X-ray selected AGN in the sample. The Alex ander et alj 
120031) X-ray catalogue is cross-correlated with the positions of the 24^m sources us- 
ing a matching radius of 3 arcsec. The X-ray emission of 44 sources in the CDF-Noith 
are likely dominated by star-formation (Georgakakis et al. 2007). These non-AGN are 
excluded from the cross-correlation. (4): fx is the fraction of X-ray detected AGN in 
the sample, i.e. fx - Nx/NjoT- 



the high redshift end of the mid-IR wedge AGN candidates. Fig. 
3 plots the total IR luminosity distribution of these sources. Most 
of the mid-IR wedge AGN candidates belong to the Luminous 
Infrared Galaxy class (LlGs), with Ltot > lO^L©. 

Mid-IR wedge AGN candidates that lie close to or are asso- 
ciated with an X-ray detection are excluded from the stacking 
analysis to avoid contamination of the stacked signal by the X- 
ray photons in the wings of the Chandra PSF The stacking re- 
sults for the remaining mid-IR wedge AGN candidates are pre- 
sented in Table |2] We estimate a soft mean X-ray spectrum for 
these sources, consistent with F x 2.1. The mean 1.5-4.0 keV 
band X-ray luminosity of this population is 6 x lO'*" erg s ' . For 
comparison, there are 670 24//m sources outside the Stern et al. 
wedge which do not lie close to an X-ray detection. Stacking 
analysis shows that the mean hardness ratio of this population is 
consistent with F a; 1 .9 and the mean X-ray luminosity is 5 x 1 0"*" 



erg s ' , i.e. similar to that obtained for the mid-IR wedge AGN 
candidates. In passing, we note that the average hardness ratio 
of the 46 mid-IR wedge AGN which are detected at X-rays is 
-0.09 + 0.01 corresponding to F ~ 1.6. 



4.3. Red/optically faint mid-IR AGN candidates 

The red/optically-faint mid-IR AGN candidates have an X-ray 
identification rate of about 19% (Table 1). Stacking the X-ray 
photons at the positions of these 83 sources that are not detected 
individually at X-rays yields a statistically significant detection 
in both the soft and the hard energy bands (table 3). The mean 
hardness ratio corresponds to a photon index of T a; 0.8 + 0.2. 
However, as the hardness ratios cannot differentiate between an 
intrinsically flat and an absorbed spectrum, it is equally possible 
that we are viewing AGN with absorbing column densities of 
~ 8 X 10"^ cm"^ for an intrinsic power-law X-ray spectrum with 
r = 1.9 and a mean redshift of z = 2. For comparison, the hard- 
ness ratio of the 20 red/optically-faint sources with X-ray detec- 
tion is 0.22+0.015 corresponding to F ~ 0.8. After the exclusion 
of one source which dominates the hard X-ray signal we obtain 
a hardness ratio of 0.09 ± 0.02 corresponding to F ~ 1.1. Only 
a small subsample of the red/optically-faint mid-IR AGN candi- 
dates, 31 out of 103, have optical counterparts. The mean photo- 
metric redshift for these 3 1 sources is Zphoio ~ 1 -6 and the mean 
IR luminosity is Ltot ~ lO'^L©, i.e. Ultra-Luminous Infrared 
Galaxies (ULIRGs). We note that one of the optically identified 
red/optically-faint mid-IR AGN candidates is a sub-mm source 
at a spectroscopic redshift of z = 2.015 (Alexander et al. 2005). 



5. Discussion 

5.1. Mid-IR wedge AGN candidates 

The X-ray bac kground population synthesis models of 
iGilh et al] (|2007|) predict a steeply rising number of Compton- 
thick AGN just below the flux limit of the CDF-N. In partic- 
ular, they predict that ~25% of the sources in the 2 - lOkeV 
band are Compton-thick AGN at fluxes ~ 10~'^er gcm~^s~' . 
Extra polating the AGN number-count distribution of Bauer et al.l 
(12004 ) yields about 70 Compton-thick AGN in the common 
Spitzer/Chandra strip at these flux levels. At fluxes of 10"'^ 
erg cm"^ s"' , comparable to the flux of the mid-IR wedge AGN, 
th is number r i ses to 200. The application of the colour criteria 
of IStern et all ( l2005h in the GOODS North sample have identi- 
fied 131 AGN with no X-ray counterpart. The number of mid-IR 
wedge AGN candidates is of the same order of magnitude as that 
predicted by the X-ray background synthesis models. We caution 
that the above is only a very coarse comparison as the XRB syn- 
thesis models cannot accurately predict the number of Compton- 
thick AGN at mid-IR wavelengths. This is because some sources 
can be so obscured that the nuclear flux will emerge only at the 
longest rest frame wavelengths, outside the IRAC bands, (Silva 
et al. 2004, Ballantyne & Papovich 2007). Moreover, at high red- 
shift iz > 2) the IRAC bands will sample only the rest frame 
galaxy emission. 

The steep mean X-ray spectral properties of these sources, 
F » 2.1, are inconsistent with the flat spectral index expected 
for Compton-thick AGN in the Chandra passband (F ~ 1; e.g. 
Georgantopoulos et al. 2007). This makes the case for a substan- 
tial fraction of Compton-thick sources among the mid-IR wedge 
AGN candidates difficult. Further clues on the nature of these 
sources can be obtained from the ratio between X-ray and in- 
frared luminosities. This is estimated Lx/^ir ~ 2x10"^, approxi- 



mately 2 dex lower than unobscured AGN (Elvis et al. 1994), and 
comparable to the ratio of Lx/^ir ~ 4 x 10"^, determined for the 
mid-IR sources outside the Stern et al. wedge (see Table 2). For 
comparison, local star- forming galaxies have L^ILm » 8 x 10"^ 
(adapted from Ranalli et al. 2003). This is a factor of few higher 
than the value derived for the mid-IR wedge candidate AGN. 
The very low value of Lx/^ir suggests that the mid-IR wedge 
candidate AGN sample is largely contaminated by normal galax- 
ies. This is demonstrated in Fig. 1, where the colour-colour track 
of the star-forming galaxy Arp 220 overlaps with the Stern et al. 
wedge at specific broad redshift intervals. 



5.2. Comparison with mid-IR power-law methods 

Mid-IR power-law selection criteria have been proposed to pro- 
duce AGN samples with a small level of galaxy contamination. 
The mid-IR power-law sources form a s ubsample of the mid- 
IR wedge sources of Stern et al. (2005). ^ Alonso-Herrero et al.l 
(2006) selected 24yum sources in the CDF-S which present red 
power-law SEDs in the Spitzer (3.6-8.0/im) bands (/v oc v" with 
a < -0.5). Nearly half of their 92 sources are not detected in X- 
rays. iDonley et al.l (l2007h apply similar power-law SED criteria 
to the Spitzer sources of CDF-N with the difference that they do 
not require their sources to be detected in 24 |/m. They consider 



47 AG N within an off-axis angle of <10 arcmin. Donl ev et al. 



('2007') find that of these, 30 are detected in the Alexander et al 



(2003) catalogue, 10 more are detected at very low significance 
(> 2.5cr) at X-ray wavelengths, and 7 remain undetected. In Fig. 
4 we show the position of the power-law selected sources on the 
mid-IR colour-colour diagram. The vast majority fall within the 
Stern wedge. We note that the 17 sources (10 low significance 
and the seven undetected) of Donley et al. generally avoid the 
central part of the CDF-N. This suggests that these have not been 
detected in the lAlexander et al](l2003.) catalogue simply because 
the exposure time is lower in these regions due to vignetting. 

In order to compare their work with our results, we estimate 
the co-added signal of the 17 sources. The results are presented 
in Table 4. We detect a statistically significant signal in both the 
soft (0.3-1.5 keV) and the hard (1.5-4 keV) band. The average 
flux in the 0.5-2 keV band is 3.4 x 10"'^ ergcm^^s"' , very 
close to the flux limit of the current CDF-N observation. The 
hardness ratio is -0.07 + 0.08 translating to a photon index of 
Fa; 1.55 ± 0.15, or alternatively a column density of Nh ~ 10^^ 
cm"^ at z=2 (for F = 1 .9). We note that when we stack only the 
7 sources which are not detected at X-rays, we obtain no sig- 
nal in either the soft or the hard band. In any case, the number of 
missed AGN in X-r ay surveys recovered with mid-IR power-law 
methods is not high. lDonlev et al.1 (|2007|) find 7 mid-IR AGN not 
detected in X-rays which form a miniscule fraction of the few 
hundreds X-ray detected AGN in the same area. This figure is 
at least one order of magnitude lower than the estimates of the 
X-ray background population synthesis models, which predict 
large numbers of Compton-thick AGN at fluxes just fainter than 
the CDF-N (Gifli et al. 2007). 



5.3. Red/optically faint mid-IR AGN candidates 

Red/optically faint mid-IR AGN candidates are more promising 
in terms of number density for being the Compton-thick pop- 
ulation predicted by the models of Gilli et al. (2007). For the 
red/optically-faint mid-IR AGN candidates we estimate a mean 
hardness ratio of 0.22 ± 0.06, which is comparable to that de- 
rived by Fiore et al. (2008; HR a 0.1) for similarly selected 



Table 2. Stacking Results for the wedge selection of Stern et al. 



Type 
(1) 



No 
(2) 



z exposure 
(3) (4) 



Soft counts 

(5) 



Soft Flux 
(6) 



1 X 10^" 
1 X 10-" 



(7) 



6x10^ 
5 X 10'"' 



Hard counts 
(8) 



Hard Flux 

(9) 



7x 10"" 
1 X 10-" 



HR 

(10) 



UUr 

(11) 

TxTcP" 
4 X 10-^ 



AGN 
Galaxies 



126 
670 



1.44 
1.05 



190 
1010 



244: 

1685 4 



: 37 (7.1a-) 
121 (15.0cr) 



147: 
1136: 



: 40 (3. So-) 
131 (8.90-) 



-0.25 : 

-0.19: 



0.05 
0.02 



The columns are: (1) sample definition; (2): Number of sources used in the stacking; (3): mean redshift; (4) effective exposure time in Ms; (5): 
net counts in the 0.3- 1.5 keV band and significance; (6): fluxin the 0.3- 1.5 keV band in units of erg cm"^ s"' ; (7): luminosity in the 0.3-1. 5keV 
band in units of ergs"' ; (8) net counts in the 1.5-4 keV band and significance; (9) flux in the 1.5-4 keV band in units of erg cm"- s"' ; (10) 
Hardness Ratio; (11): Soft X-ray (1.5-4 keV) to IR luminosity ratio 



Table 3. Stacking results for the red/optically faint AGN candidates 



Type 
(1) 



No 
(2) 



exposure 

(3) 



Soft counts 

(4) 



Soft Flux 

(5) 



Hard counts 
(6) 



Hard Flux 

(7) 



HR 

(8) 



AGN 83 



125 



105±29 (3. So-) 5.5 x 10" 



164±33(5.4cr) 1.3 x 10" 



0.22 + 0.06 



The columns are: (1) sample definition; (2): Number of sources used in the stacking; (3)) effective exposure time in Ms; (4): net counts in 
the 0.3-1.5 keV band and significance; (5): flux in the 0.3-1.5 keV band in units of ergcm"^ s"' ; (6) net counts in the 1.5-4 keV band and 
significance; (8) Hardness Ratio. 



Table 4. Stacking results for the mid-IR power-law methods 



Type 
(1) 



No 
(2) 



exposure 

(3) 



Soft counts 

(4) 



Soft Flux 

(5) 



Hard counts 
(6) 



Hard Flux 

(7) 



HR 

(8) 



AGN 17 



18 



78±13(7.8o-) 3.1x10" 



68±14 (5.90-) 3.7 x 10" 



-0.07 ± 0.08 



The columns are: (1) sample definition; (2): Number of sources used in the stacking; (3)) effective exposure time in Ms; (4): net counts in 
the 0.3-1.5 keV band and significance; (5): flux in the 0.3-1.5 keV band in units of ergcm"^ s"' ; (6) net counts in the 1.5-4 keV band and 



significance; (8) Hardness Ratio. 
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Fig. 4. The mid-IR power-law selected sources of Donley et al. 
(2007); filled symbols are the X-ray detections, open circles the 
X-ray non detections and crosses the marginal 2.5 sigma detec- 
tions (see §7.2). 



sources in the CDF-South. Our hardness ratio corresponds to a 



photon index of T « 0.8 ± 0.2, in the range of the spectral index 
of Compton-thick AGN which are dominated by the reflected 
component. The mid-IR-to-X-ray flux ratio of these sources is 
F(24//m)/F(2 - 8keV) ~ 6 x 10'*'[//Jy]/[erg cm"^ s"'], sim- 
ilar to that of the prototype Compton- thick AGN NGC 1068 
at z a! 2 (Daddi et al. 2007). It is therefore likely that a large 
fraction of the red/optically-faint AGN candidates are associ- 
ated with Compton-thick sources at high redshift. Alternatively, 
they may be low-luminosity AGN obscured by column densi- 
ties, Nh ~ 8 X 10^^ cm"^ at redshift of z=2. Good quality X-ray 
spectra could diflFerentiate between the above scenaria. However, 
the average flux of these sources is about a factor of four lower 
than the current flux limits of the 2 Ms CDF-N in both the 0.3-1 .5 
and 1.5-4 ke V bands. This means that a total of 8 Ms of Chandra 
observations are needed just to detect this population in the soft- 
band at least in the central part of the field-of view, where the 
background will remain photon limited (see Alexander et al. 
2003). Alternatively, near-IR spectra could be useful in detecting 
AGN signatures such as a strong [OIII] Line. 

The mid-IR colours of the red/optically-faint mid-IR AGN 
candidates provide additional information on their nature. In Fig. 
1 the majority of these sources lie in the region of the param- 
eter space occupied by ULIRGs, such as Arp220 or Mrk273, 
at redshifts z » 2, i.e. galaxies in formation, possibly through 
mergers (Mirabel & Sanders 1996). This is consistent with the 
mean photometric redshift of » 1 .6 and the mean total infrared 
luminosity of k, lO'^L© estimated for the 31 red/optically-faint 
mid-IR AGN candidates with optical counterparts. High reso- 
lution HST/ACS images also suggest that these sources include 
interacting systems. There are 20 sources bright enough in the ; 



or z bands for visual inspection. Many of them (13) either show 
disturbed optical morphology or are in close pairs/associations. 
One of them is associated with a sub-mm source at z = 2.015. 
Although the morphological subsample represents only about 
20% of the red/optically-faint mid-IR AGN candidates, the evi- 
dence above is suggestive that this population includes galaxies 
captured in the process of formation via tidal interactions. In this 
respect, these obscured AGN bear many similarities to high red- 
shift systems believed to undergo the concurrent growth of their 
stellar population and their central SBH. Both sub-mm galaxies 
at z « 2 (Alexander et al. 2005) and mid-IR excess 24yum se- 
lected sources at 1.4 < z < 2.5 (Daddi et al. 2007) are suggested 
to harbour Compton-thick AGN, in addition to star-formation 
activity. 

Interestingly at lower redshift, z ^ 1, most of the ob- 
scured AGN are associated with early-type galaxies rather than 
systems in formation (e.g. Rovilos & Georgantopoulos 2007). 
Georgakakis et al. (2008) used stacking analysis to search of 
heavily obscured AGN among different types of galaxies in the 
range 0.4 < z < 0.9. They showed the emergence of an X-ray 
faint and obscured AGN population, possibly including some 
Compton-thick sources, among galaxies in the red cloud, where 
most of the early-type quiescent systems are found. The evi- 
dence above may suggest a change in the properties of the galax- 
ies that host obscured, possibly Compton-thick, AGN. At high 
redshift, z ~ 2, an increasing body of evidence points to a link 
with galaxies in formation, while at lower redshifts, z ^ 1, they 
are found in early type hosts. 

6. Conclusions 

We investigate the X-ray properties of mid-IR selected AGN in 
the CDF-N, using the Spitzer 24fj.m sample. Our conclusions can 
be summarised as follows. 

- Stacking analysis of the mid-IR wedge AGN candidates re- 
veals a soft mean X-ray spectrum, consistent with T x 2.1, 
and a low X-ray-to-IR luminosity ratio, both suggesting that 
the sample is contaminated by normal galaxies. We conclude 
that the mid-IR colour 'wedge' methods fail to easily iden- 
tify AGN, when applied to faint optical magnitudes. 

- Mid-IR 'power-law' selection techniques are more success- 
ful in finding AGN. These sources overlap substantially with 
X-ray selected AGN however, and do not represent a new X- 
ray weak population. Moreover, the total number of mid-IR 
'power-law' AGN in the CDF-N is well below the total num- 
ber of Compton-thick AGN predicted by X-ray background 
population synthesis models. 

- Sources with high 24|^m-to-optical {ful fopt > 1000) flux 
ratio and red colours {Rab - ^3.6 > 3.7) have hard stacked X- 
ray spectrum (F « 0.8), suggesting Compton-thick systems. 
Low Luminosity AGN with column densities of 8 x 10^^ 
cm"^ or a combination of the two. Most of these sources 
occupy a distinct region of the colour-colour mid-IR dia- 
gram and lie outside of the Stern wedge. The mid-infrared 
colours and luminosities of these sources are consistent with 
ULIRGs at z ~ 2. For the optically bright subsample of this 
population, the HST/ACS images show evidence for interac- 
tions, suggesting systems in the process of formation. 
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